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Diffusion kinetics in contact of electromelted corundum refractories with the glass melt was investigated. The 
diffusion coefficients of Al 3+ and their temperature dependence were determined. The basic tenets of diffusion 
kinetics in service of refractories in glass furnaces and feeder channels are examined. 


Electromelted refractories are predominantly used in the 
glass industry. Corundum refractory articles are highly re¬ 
commended in lining for the feeder channels of glass mold¬ 
ing machines. Their high stability at the temperatures of 
preparation of the glass melt for molding ensures a long life¬ 
time and prevents the appearance of glass flaws in the fin¬ 
ished products. 

At the same time, the complex process of the reaction of 
refractories with the glass melt, including surface phenom¬ 
ena, chemical reactions, transfer phenomena (diffusion, vis¬ 
cous flow, heat transfer) and other factors has been insuffi¬ 
ciently investigated. This especially concerns diffusion phe¬ 
nomena; the information on the diffusion coefficients of dif¬ 
ferent ions is limited and primarily concerns diffusion of Na + 
ions at low temperatures [1]. 

Experimental methods. Samples of electromelted co¬ 
rundum refractories were used for the study (mass content, 
%): 1) 100 A1 2 0 3 ; 2) 96.50 A1 2 0 3 , 2.65 Si0 2 , 0.41 Na 2 0, 
0.05 K 2 0, 0.12 Fe 2 0 3 , 0.27 other oxides. Composition of the 
glass (%): 75 Si0 2 , 10 CaO, 15 Na 2 0. The experiments con¬ 
sisted of melting glasses of the given composition, obtaining 
data on the viscosity and density of the glasses with consid¬ 
eration of dissolution of the refractories, determination of the 
concentration ratios in dissolution, the diffusion characteris¬ 
tics, etc. 

To study corrosion of the refractories in the glass melt, 
different static and dynamic methods were used [2] which al¬ 
lowed obtaining a comprehensive evaluation of the glass sta¬ 
bility of the refractory materials without separating the indi¬ 
vidual constituents of the process. The methods are relatively 
graphic in conducting comparative tests of samples of 
refractories using the same experimental setup. Comparable 
results cannot be obtained on different setups, so that theore¬ 
tically substantiated methods of diffusion kinetics in moving 
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media for different conditions were used to study the diffu¬ 
sion characteristics in reaction of the refractories with the 
glass melt. 

Such methods are the rotating disk method (sample of 
the investigated refractory) ensuring the perpendicularity of 
diffusion flow to the surface [3, 4]. When this method is 
used, the following picture of movement of the liquid is ob¬ 
served. Ascending flow to the disk is observed at a sufficient 
distance from the rotating disk, and in the layers immediately 
adjacent to the disk, the liquid acquires rotational motion, 
and the angular velocity increases in approaching the disk, to 
a value equal to the angular velocity of the disk itself. The ro¬ 
tating disk has an important feature that ensures equal acces¬ 
sibility of the surface with respect to diffusion, which is very 
important for real solids whose surface is usually not homo¬ 
geneous. In addition, the thickness of the boundary layers in 
this case will have a constant value over the entire surface. 

The analysis of the hydrodynamics of movement of a liq¬ 
uid given in [3] allowed obtaining an equation for calculating 
the rate of dissolution of solid samples with consideration of 
the diffusion coefficients: 

j = 0.62 D 2//3 v “ (c s - c) co 1 / 2 , 

where D is the diffusion coefficient; v is the kinematic vis¬ 
cosity; (c s - c) is the difference in the saturation concentra¬ 
tions at a sufficient distance from the rotating sample; co is 
the angular velocity of rotation of the sample. 

The high-temperature versions of the rotating sample 
method were used by P. M. Shurygin and other scientists in 
investigating diffusion kinetics for metallurgical melts (slags, 
salts) [5, 6]. The research of M. Trukhlarova and O. Veprek 
(1969 - 1970) should also be noted, as they used the rotating 
quartz cylinder method in studying diffusion of Si0 2 in so¬ 
dium-calcium-silicate glass [1]. B. Ya. Ten [7] used the 
method of calculating the diffusion coefficient of Si0 2 with 
the kinetic curves of dissolution of silica grains. 
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In conducting our experiments, the samples were at¬ 
tached in the corresponding holder and rotated with low an¬ 
gular velocities, which excluded turbulent modes in consid¬ 
eration of the high viscosity of the melts. The diffusion coef¬ 
ficients were calculated with the experimental data using 
the above equation. The saturation concentration (solubility 
of A1 2 0 3 in the glass melt) was determined with the full 
Na 2 0 - CaO - A1 2 0 3 - Si0 2 diagram. 

Chemical reaction of refractories with the glass melt. 
The reaction of a glass melt with a refractory, including the 
diffusion kinetics, is complex and requires special studies 
with analysis of an important amount of experimental data. 
Nevertheless, based on the available data, it is possible to hy¬ 
pothesize a two-stage model of the reaction with consider¬ 
ation of structural factors and diffusion kinetics [4, 8]. 

On contact of the glass melt with the refractories, the re¬ 
fractory first reacts with the most active alkali ions in the 
glass melt. The scheme of the reaction can be represented as 
follows. On “attack” of the alkali ions on the crystal lattice of 
the refractory material, it is destroyed with subsequent mi¬ 
gration of the constituent elements into the glass melt. In the 
case of corundum refractories, Al 3+ ions and oxygen are 
formed, and for oxygen, bridge ions are converted into free 
nonbridge ions: 

20 ° -> 0 2 “ + 0 “ 

These ions diffuse into the glass melt and participate in 
formation of coordination polyhedrons. It is necessary to 
keep in mind here coupling of the flows of charged particles. 
Departing from conditions of electroneutrality, movement of 
ions of the same sign should be accompanied by movement 
of ions of the opposite sign. In this case, the diffusion kinet¬ 
ics will change as a function of the chemical and structural 
changes in the glass in the contact zone. 

In the initial period of the reaction, the driving force of 
diffusion of the Al 3+ ion is much greater than for oxygen ions 
due to the concentration ratios; on the contact boundary, the 
concentration of Al 3+ is maximum and equal to zero in the 
glass, while the difference in the concentrations of oxygen 
ions is insignificant (in the glass, the concentration of active 
oxygen ions is determined by the Na 2 0 and CaO content). In 
consideration of the size of the Al 3+ (0.57 A) and oxygen 
ions (1.36 A), the driving force of diffusion of Al 3+ will be 
much higher than for oxygen ions. 

In migration of the Al 3+ ion out of a corundum refractory 
into glass, its coordination number for oxygen changes from 
6 to 4. In glass, the Al 3+ ion, surrounded by four oxygen ions, 
can isomorphically replace the Si 4+ ion in silicon-oxygen 
formations. Highly coherent single silicon-aluminum-oxy¬ 
gen formations arise. This is indicated by the increase in the 
viscosity and surface (interfacial) tension. As indicated in 
[8], the state of structural oxygen in the glass-forming melt is 
determined by the equilibrium conditions: 

0 ° + 0 2 “ 20 “ 


These conditions will change when the 0 2 “ and 0“ ions 
arising in decomposition of the crystal lattice of corundum 
appear. Bridge 0° ions in the glass-forming melt are already 
bound with silicon and do not participate in formation of 
ionic aluminum-oxygen polyhedrons (tetrahedrons), so that 
the excess of active oxygen ions, predominantly 0“, is used 
to construct these tetrahedrons. The use of the oxygen ions in 
the glass for this purpose will approximately continue until 
the condition (Na 2 0 + Ca0)/Al 2 0 3 = 1 is satisfied. Further¬ 
more, oxygen ions formed in decomposition of the crystal 
lattice of the corundum refractory will be used for building 
aluminum-oxygen tetrahedrons until the saturation concen¬ 
tration of A1 2 0 3 is attained in the glass. The driving force of 
diffusion for oxygen ions should increase due to the impor¬ 
tant difference in the concentrations of active oxygen ions on 
the contact boundary and in the glass and should decrease for 
Al 3+ because of a decrease in the difference in concentra¬ 
tions. However, the difference in the size of the ions can 
level out the diffusion characteristics of the ions. The overall 
rate of the diffusion process will slow, however. 

In the contact zone, saturation of the glass melt with 
A1 2 0 3 occurs according to the Na 2 0 - CaO - A1 2 0 3 - Si0 2 
phase diagram. Applying the composition of the glass used 
in the diagram and joining it with the A1 2 0 3 apices, we ob¬ 
tain the direction of the reaction of the corundum refractory 
with the glass melt and the compositions of intermediate and 
crystallized products. As the diagram suggests, these prod¬ 
ucts are mixtures corresponding to minerals of the sodium 
and calcium feldspar series. Products of the sodium feldspar 
series are primarily formed due to the higher reactivity of the 
Na + ions and in particular, nepheline and carnegieite: 

2A1 2 0 3 + Na 2 0 + 2Si0 9 = Na 2 0 • A1 2 0 3 • 2Si0 2 + A1 2 0 3 . 

nepheline, carnegieite 

A multitude of small crystals of these compounds is 
founding the contact zone in the glass, which has also been 
noted by other investigators [2]. 

Features of diffusion kinetics. High-temperature reac¬ 
tions of a highly viscous melt with a solid and correspond¬ 
ingly diffusion phenomena are very complex. 

Transfer of a substance in a moving liquid is due to two 
different mechanisms. First, when there is a concentration 
difference in the liquid, molecular or ionic diffusion arises; 
second, particles of the substance dissolved in the liquid are 
entrained by the liquid as it moves and are conveyed along 
with it. The special processes characterize convective trans¬ 
fer of matter in a liquid. 

It should be noted that the substantiated theoretical as¬ 
sumptions of convective transfer were basically elaborated 
for low-viscosity liquids, and the necessity of studying diffu¬ 
sion kinetics in highly viscous liquids was emphasized in 
[3,4], since the experimental data for them are clearly insuf¬ 
ficient. 

The Reynolds number and diffusion number, which re¬ 
spectively characterize the movement mode and the physico- 
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Fig. 1 . Character of the change in the concentration of A1 2 0 3 (7) 
and rate of melt movement (2) in a rotating disk sample in the direc¬ 
tion of diffusion flow: v ) rate of melt movement on the inner boun- 
dary of the hydrodynamic layer; v) same for the outer boundary. 


chemical properties of the liquid, are the basic parameters of 
a moving liquid. Movement of a glass melt in glass furnaces 
is characterized by low values of the Reynolds number and a 
regime of developed turbulence is difficult to create in ge¬ 
neral for a highly viscous glass melt. For this reason, the im¬ 
portance of the Prandtl diffusion number increases as one of 
the basic characteristics of a highly viscous melt: 

Pr = v/D. 

For glass melts, the Prandtl diffusion number can be 
equal to 10 8 - 10 10 m 2 /sec. Such high values emphasize some 
features of the diffusion process in highly viscous melts. 

For experimental and industrial conditions, the entire 
volume of a highly viscous melt can arbitrarily be divided 
into three regions: 

the region of the basic volume of the melt in which the 
concentration of dissolved refractory matter is insignificant 
in comparison to the saturation concentration; movement of 
the glass melt in this volume is due to organized convective 
diffusion (experimental vessel) or production and convective 
flows (glass furnace); 

the region of constant concentration distribution and con¬ 
vective diffusion at a certain distance from the surface of the 
refractory - the hydrodynamic boundary layer (in the scien¬ 
tific literature, it is also called the “Prandtl layer”); laminar 
movement occurs here and the viscosity of the melt plays the 
basic role; in studying the diffusion characteristics by the ro¬ 
tating disk method and low values of the Reynolds number, 
the thickness of the hydrodynamic layer is commensurate 
with the radius of the disk [3]; it was noted in [9] that in dis¬ 
solution of refractories in melted glass and other silicate 
melts that are highly viscous, not very mobile liquids, rela¬ 
tively thick boundaries layers which can attain 10“ 2 m and 
more are formed; 

the region of a rapid change in concentration and diffu¬ 
sion of directly near the surface of the refractory material - 
the diffusion layer (also called the “Nemst layer” in the sci¬ 
entific literature); saturation concentrations are attained here; 
for highly viscous glass melts enriched with high-melting 
components of the refractories, this layer can be considered 
strongly bound with the surface of the refractories; the prod¬ 


ucts of dissolution migrate into the bulk of the melt by con¬ 
vective diffusion through the hydrodynamic boundary layer; 
the products of dissolution go from here into the bulk of the 
melt by convective diffusion through the hydrodynamic 
boundary layer; the thickness of the diffusion layer 5^ is 
small and in agreement with the data in [3,4], can be approx¬ 
imately determined from the relation: 

8 0 = 8/Pr‘/ 3 , 

where 8 is the thickness of the hydrodynamic layer. 

It follows from this relation that the thickness of the dif¬ 
fusion layer decreases significantly with an increase in the 
viscosity of the liquid (for glass melts, at a hydrodynamic 
layer thickness of 10“ 2 m and a Prandtl number of 10 9 , the 
thickness of the diffusion layer will be 10 _5 m, while the 
Prandtl number will only be 10 3 and correspondingly, the dif¬ 
fusion layer thickness will be 10 “ 3 m for aqueous medium 
with the same hydrodynamic layer thickness). Similar results 
were obtained in [10] for the reaction of a glass melt with 
refractories (with consideration of the assumptions made and 
experimental conditions). 

The generalized dependences of the change in the re¬ 
ported concentration of A1 2 0 3 and the rate of movement of 
the melt in the direction of diffusion flow perpendicular to 
the surface of the rotating samples (axis v) are shown in 
Fig. 1. The thickness of the diffusion and hydrodynamic lay¬ 
ers is determined by the points of intersection of tangents to 
the curves with the level of limiting values of the concentra¬ 
tion and movement rate. These values are attained on the in¬ 
ner and outer boundaries of the hydrodynamic layer. As 
noted, beyond the limits of the outer boundary of the hydro- 
dynamic layer, the character of movement of the liquid 
changes, related to the presence of ascending flow to the 
disk. 

The experimental dependences of the rate of dissolution 
of the disk samples on the angular velocity of rotation and 
the temperature are shown in Fig. 2. The linear character of 
the dependence indicates a diffusion mode of dissolution of 
the refractory samples, i.e., the dissolution rate is not limited 
by the rate of the chemical reaction, but by removal of the 
products formed from the reaction zone. The temperature de¬ 
pendence of the diffusion coefficient is shown in Fig. 3. The 
activation energy of the diffusion process was determined 
with the experimental data obtained and was 258 kJ/mole. 
The low values of the diffusion coefficients and high activa¬ 
tion energy of the process indicate the elevated resistance of 
electromelted corundum refractories to sodium-calcium-sili¬ 
cate glass melts. 

Simlar experiments were conducted with samples of in¬ 
dustrial corundum refractories containing small amounts of 
silica, sodium, potassium oxides, etc., in addition to A1 2 0 3 . 
The content of these oxides increases the thermal stability of 
the refractories, but their glass stability decreases relatively 
significantly. This is indicated by the values of the diffusion 
coefficient of Al 3+ , whose value increases by more than one 
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Fig. 2. Rate of dissolution of refractory sample (A1 2 0 3 ) on angular 
velocity of rotation at temperatures of 1350°C (7), 1400°C (2), 
1450°C (3), and 1500°C ( 4 ). 

order of magnitude, especially at high temperatures. This is 
due to the presence of a small amount of glass phase contain¬ 
ing alkali oxides. Penetration of Na + from the glass into the 
refractory also takes place over the glass phase. 

The increased mobility and reactivity of the sodium ions, 
especially at high temperatures, create the conditions for de¬ 
struction of the crystal lattice of the refractory from the in¬ 
side and significantly increase the diffusion coefficient of 
Al 3+ . For this reason, it is best to use electromelted corundum 
refractories for the feeder pans of glass-forming machines, 
where the temperature is 150 - 200°C lower than in the melt¬ 
ing part of the glass furnace. The comparatively low temper¬ 
ature in the feeders and the related high viscosity of the glass 
melt decrease the aggressiveness of the melt and the diffu¬ 
sion coefficient of the ions. At the same time, the increase in 
the rate of movement of the glass melt in the feeder channels 
causes accelerated removal of the products formed from the 
zone of the chemical reaction. 

Diffusion kinetics in industrial conditions. In consider¬ 
ation of the importance of the corrosion resistance of refrac¬ 
tory materials in glass furnaces, the reaction of these materi¬ 
als with the glass melt has been investigated in many studies, 
most of which have been published [2, 11]. 

Let us consider some features of diffusion phenomena on 
contact of the glass melt with the walls of feeder channel 
pans made of electromelted corundum refractory. Similar 
phenomena will also take place on contact of the glass melt 
with the wall blocks of the glass furnace (the character of the 
diffusion kinetics with consideration of the Zr0 2 content is 
not examined here). 

A diagram of movement of a melt and distribution of the 
rates in flow of a glass melt moving relative to the immobile 
refractory wall parallel to axis y (production flow) is shown 


-log D 



Fig. 3. Temperature dependence of the diffusion coefficient in dis¬ 
solution of A1 9 0 3 in a glass melt. 




Fig. 4. Character of diffusion kinetics in an immobile refractory: 
a) diagram of movement of the melt and velocity distribution in flow 
of a glass melt in convective diffusion; b ) character of the change in 
the concentration of A1 2 0 3 (7) and rate of melt movement (2) in the 
direction of diffusion flow. 


in Fig. 4 a. The flow rate increases from zero (on the surface 
of the refractory) to defined values within the limits of the 
hydrodynamic layer. Diffusion takes place in the direction 
perpendicular to the direction of movement of flow of the 
glass melt, i.e., along axis x. The basic process of diffusion of 
ions takes place in the thin diffusion layer and beyond its 
boundaries, convective diffusion takes place in the hydrody¬ 
namic layer. 

The generalized dependences of the change in the pro¬ 
duction concentration and A1 2 0 3 and the rate of movement of 
the glass melt are shown in Fig. 4b. In this case, in contrast to 
the rotating samples, the rate on the surface of the immobile 
refractory (tank block, feeder walls) is equal to zero. 

It should be noted that the phenomena examined are sec¬ 
ondary. Formation of a diffusion layer primarily takes place 
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due to the chemical reaction, concentration gradient, and sur¬ 
face forces, and the character of the change in them is given 
in [11, 12]. 

We should also note that in the refractory-glass melt 
contact zone, the portions of glass melt enriched with A1 2 0 3 
and poor in Na 2 0 drop down and new portions of glass melt 
with a lower A1 2 0 3 content or with no A1 9 0 3 appear in their 
place (gravitational convection). The surface forces acting 
within the boundaries of the diffusion layer significantly af¬ 
fect dissolution of refractories, despite their low intensity. 
These forces can change the thickness of the diffusion layer. 
When it decreases, entrainment of dissolved matter (A1 2 0 3 ) 
from the reaction zone is accelerated due to convective diffu¬ 
sion, corrosion of refractories is intensified, and flaws appear 
in the glass melt in the form of ripples, stones, etc. 

A difference in the size of surface (interphase) strains on 
the boundary of the diffusion and hydrodynamic layers. A 
smaller difference decreases the rate of movement of the 
glass melt in the contact zone and its effect on the diffusion 
layer by decreasing the corrosion rate of the refractory. This 
is especially marked in corrosion on the surface of the glass 
melt and in vertical pit corrosion [11]. 

The effect of surface and concentration factors in use of 
chromium oxide refractories in melting green glass colored 
with Cr 2 0 3 (chromium oxides are surface-active components 
in the glass) is marked. In this case, the difference in the sur¬ 
face (interphase) stresses is insignificant, and the saturation 
concentration and concentration gradient (c s - c ) for Cr 2 0 3 
will be much smaller than for A1 2 0 3 . This significantly de¬ 
creases the corrosion rate of the refractories. The decrease in 
the diffusion coefficient also causes an increase in the viscos¬ 
ity of the glass melt in the contact zone, which is usually ob¬ 
tained by cooling the corresponding section of the refractory 
lining, for example, the top row of wall blocks at the level of 
the surface of the glass melt and the neck blocks. When 
Bakor refractories are used, the viscosity increases due to 
combined dissolution of Zr0 2 and A1 2 0 3 . 

The basic factors of diffusion kinetics in the reaction of 
corundum refractories with the glass melt, including in in¬ 
dustrial conditions, were thus revealed. The experimental re¬ 


sults obtained indicate the diffusion character of dissolution 
of corundum refractories in a glass melt. The diffusion coef¬ 
ficients of Al 3+ ions determined with the experimental data 
are 10“ 11 - 10“ 12 m 2 /sec. 

Analysis and consideration of diffusion kinetics factors 
will increase the efficiency of use of refractories in glass fur¬ 
naces. 
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